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Abstract 

First laws of black hole mechanics, or thermodynamics, come in a variety of different forms. In 
this paper, from a purely post-Newtonian (FN) analysis, we obtain a first law for binary systems 
of point masses moving along an exactly circular orbit. Our calculation is valid through 3FN order 
and includes, in addition, the contributions of logarithmic terms at 4FN and 5FN orders. This first 
law of binary point-particle mechanics is then derived from first principles in general relativity, and 
analogies are drawn with the single and binary black hole cases. Some consequences of the first 
law are explored for FN spacetimes. As one such consequence, a simple relation between the FN 
binding energy of the binary system and Detweiler's redshift observable is established. Through it, 
we are able to determine with high precision the numerical values of some previously unknown high 
order FN coefficients in the circular-orbit binding energy. Finally, we propose new gauge invariant 
notions for the energy and angular momentum of a particle in a binary system. 
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I. INTRODUCTION 



A. Motivation 

Inspiralling and coalescing binary systems composed of black holes and/or neutron stars 
are among the most promising sources of gravitational radiation that might be detectable by 
forthcoming ground-based interferometers, such as Advanced LIGO and Advanced Virgo, as 
well as by future space-based observatories [1]. The detection and analysis of these signals 
require very accurate theoretical predictions, for use in the construction of gravitational- wave 
templates. Given the complexity of the Einstein field equations, continued progress requires 
that we must (i) cleverly identify which parts of the problem naturally lend themselves to 
simplification, and (ii) rely on a combination of approximation and numerical methods. 

There are two main approximation schemes for studying the relativistic dynamics of 
compact binary systems, and the associated emission of gravitational radiation: the post- 
Newtonian (PN) approximation, which is well-suited to describe the inspiralling phase of 
arbitrary mass ratio compact binaries in the small velocity and weak field regime {v ^ c) [2j , 
and black hole perturbation theory, which provides an accurate description of extreme mass 
ratio binaries (mi ^ ^2), even in the strong field regime [3j. 

Each of these approximation methods frequently relies on a simplified description of one 
or both of the compact objects in terms of structureless point particles, characterized solely 
by their masses mi and m2, and eventually their spins. Even though the notion of a point 
mass has been shown to be ill-defined in the exact theory of general relativity [HIS],^ it can be 
made sensible in approximation methods such as PN expansions or black hole perturbation 
theory. This idealization is particularly convenient in order to carry out calculations up to 
very high orders, as required for gravitational-wave searches relying on matched filtering. 

Except for the occurrence of a gradual inspiral driven by radiation-reaction, the orbits 
of stellar mass compact binaries can be considered to be circular, to a very high degree 
of approximation. Mathematically, the approximation of an exactly closed circular orbit 
translates into the existence of a helical Killing vector (HKV), along the orbits of which the 
spacetime geometry is invariant. This HKV field (say K°'da = dt + where Q is the 

constant circular orbit frequency) can be viewed as the generator of time translations in a 
co-rotating frame. Given the high astrophysical relevance of this approximation, helically 
symmetric spacetimes have been studied extensively in the literature [614T6] ; these numerous 
works range from formal mathematical analyses to practical calculations of initial data for 
binary black holes and binary neutron stars. 

In general relativity, it is known that helically symmetric spacetimes cannot be asymp- 
totically flat [12]. This can easily be understood from a heuristic point of view: in order to 
maintain the binary on a fixed circular orbit, the energy radiated in the form of gravitational 
waves needs to be compensated by an equal amount of incoming radiation. Far away from 
the source, the resulting system of standing waves rapidly dominates the energy content of 
the spacetime, such that the falloff conditions necessary to ensure asymptotic flatness can- 
not be satisfied. Nevertheless, asymptotic flatness can be recovered if, loosely speaking, the 
gravitational radiation can be "turned off". There are two well-known approximations to 
general relativity for which this can be achieved: the Isenberg, Wilson and Mathews approx- 



Physically, if one tries to compress an extended body down to a single point, a black hole will form before 
the point-particle limit is reached. 
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imation (or conformal flatness condition) [T7HT9] . and the post-Newtonian approximation. 
In the latter, it is in principle possible to unambiguously disentangle the conservative part 
of the orbital dynamics, from the dissipative effects related to radiation-reaction (at least 
up to 3.5PN order^). Consequently, in the PN approximation, one can justiflably consider 
non-radiative, helically symmetric, asymptotically flat spacetimes. 



B. Overview 



In this paper we shall consider such non-radiative, helically symmetric, asymptotically 
flat spacetimes, in which both compact objects are modelled as point masses mi and m2 
moving on exact circular orbits. Building on previous works [20ll2Tj .^ we compute Detweiler's 
redshift observables zi and Z2 [22], which represent the redshift of light rays emitted from 
the particles and received on the helical symmetry axis perpendicular to the orbital plane. 
On the other hand we compute also the total ADM mass M and total angular momentum 
J of the system of two point particles. Our calculations are carried through high post- 
Newtonian order, being complete up to 3PN order and including also the leading-order 4PN 
and next-to-leading 5PN logarithmic contributions. 

The ADM mass M, angular momentum J, and redshifts Zi 2 are all functions of the 
three independant variables of the problem, namely the orbital frequency fl that is imposed 
by the existence of the HKV, and the individual masses mi and m2 of the particles. We 
shall prove from our high-order post-Newtonian calculation that the variations of the ADM 
quantities are linked to the variations of the individual masses by the flrst law of the binary 
point-particle mechanics (or "thermodynamics") 

6M - n5J = zi5mi + Z25m2. (1.1) 

We shall demonstrate that this law is actually a particular case, valid when one assumes 
the existence of a HKV covering the entire spacetime, of the generalized law of black hole 
mechanics obtained by Friedman, Uryii, and Shibata [11]. Various consequences of this 
law will also be investigated and discussed in the framework of post-Newtonian spacetimes, 
notably the interesting relation 

M - 2nj = miZi + m2Z2 ■ (1.2) 



As an application of the flrst law (1.1) we shall be able to determine the numerical values 
of some previously unknown post-Newtonian coefficients (at 4PN, 5PN and 6PN orders) in 
the binding energy of the compact binary on a circular orbit. We shall also propose some 
gauge invariant definitions for the energy and angular momentum of a point particle in a 
binary system. 

This paper is organized as follows. In Sec. [TTjwe complete the PN calculations performed in 
Papers I and II by taking into account new logarithmic contributions coming from memory- 
type hereditary terms in the metric, and show that including these new contributions yields 



the first law (1.1), which is thereby proved up to 3PN order plus the log-terms at 4PN and 



5PN orders. In Sec. Ill we show that this law is a particular case of the known general 



^ As usual we refer to nPN as the order corresponding to terms 0{c in the equations of motion. 
^ The papers [20] and [21] are hereafter referred to as Papers I and II, respectively. 
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law of mechanics for systems of black holes and extended fluid bodies. Analogies with the 
single and binary black hole cases are drawn, and various discussions and alternative proofs 
are presented. Finally, in Sec. IV these results are applied (i) to provide gauge invariant 
candidates for the perturbed energy and angular momentum of a point particle in circular 
orbit about a Schwarzschild black hole, and (ii) to the numerical determination of high- 
order PN coefficients in the circular-orbit binding energy. We end up with a list of potential 
applications which are left for future work. 



II. POST-NEWTONIAN DERIVATION OF THE FIRST LAW 

In Papers I and II, motivated by high-accuracy comparison between the post-Newtonian 
approximation and self-force (SF) perturbative results [22], we derived the redshift observ- 
ables z\ and Z2 of a circular-orbit binary up to order 3PN, and included the specific contribu- 
tions of logarithms at 4PN and 5PN orders (see also Ref. p3] for the 4PN logarithm). In the 
limit of small mass ratio, we showed a very good agreement with numerical SF computations 
based on first-order black hole perturbation theory. 



A. Logarithmic contributions coming from memory-type hereditary terms 

We have found that the previous calculation of the 4PN and 5PN logarithmic terms had 
overlooked a subtle point related to the assumption of helical symmetry: in the Appendix of 
Paper II, it was shown that a class of logarithmic terms contributing to the metric appeared 
in the form of an infinitesimal gauge transformation. These terms were left out of the 
calculation, because they could not contribute to the coordinate invariant relations computed 
for quasi- circular orbits. That conclusion is correct in the case of the physical problem 
where both conservative and dissipative effects are included, such that the binary decays 
and coalesces within a finite amount of time (from the point of view of a distant inertial 
observer). However, once the helical symmetry is imposed, the binary must be seen as 
having orbited for an infinite amount of time, and the infinitesimal gauge transformation 
used in the Appendix of Paper II turns out to become meaningless. In the present Section 
we correct for this effect and consider only a class of gauge transformations allowed by the 
HKV. We shall show that this adds extra logarithmic contributions at 4PN and 5PN orders, 
but that these do not affect the self-force regime investigated in Papers I and II. 

The logarithms at 4PN and 5PN orders are produced by so-called "hereditary" terms pi] 
when the helical symmetry is imposed. They can be computed by using the "instantaneous" 
propagator given by 




This propagator depends explicitly on the orbital frequency Vt imposed by the helical sym- 
metry, through the length scale A = 27rc/i7. Indeed a certain Finite Part (FP) operation 
at -B = is required (see e.g. Papers I and II), which involves the regulator {r/X)^ and 
yields when i? — )■ some logarithms of the type ln(r/A). At quadratic non-linear order these 
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logarithms are contained into the following piece of the "gothic" metric:^ 



5hf 



J- 



1 [AM (2) k'^k^ 
V c c 



(2.2) 



The first term in the right-hand-side (RHS) corresponds to gravitational-wave tails (back- 
scattering of linear waves onto the background curvature associated with the monopole M 
of the source), while the second term is due to the non-linear memory effect (radiation of 
waves by the stress-energy distribution of linear waves). The tensor z"^(n, u) in Eq. (2.2) 



represents the coefficient of the dominant 1/r term of the (non-static part of the) linearized 
metric h^^ at future null infinity, i.e. when r — +oo with the retarded time u = t — r/c 
kept fixed. It is in the form of a sum of multipoles parametrized by mass-type moments 
Ml{u) and current-type moments Sl{u): 



^00 



1^2 



1^2 



e.^2 



C^+3(£+l)! 
2i 



u 



(2.3a) 
(2.3b) 
(2.3c) 



Our notation for multi-indices such as L = ii ■ ■ ■ ii (and so on) is the same as in Paper II. 
The superscript {£) refers to time derivatives. In the second term of the RHS of Eq. (2.2), 
= (1,?^) is the outgoing Minkowskian null vector and a{n,u), which is essentially the 
energy carried away by the waves in the direction n = x/r, is given by 



(2.4) 



The whole computation of the log-terms in Paper II was based on the first contribution 



(tail) in Eq. (2.2). The second contribution (memory) was discarded because the log-terms 
therein are formally in the form of an infinitesimal gauge transformation. However, the argu- 
ment overlooks the fact that the needed gauge transformation involves some anti-derivatives 



of the quantity a defined by (2.4); of. the Appendix of Paper II. Since the average of a 



over all directions represents the energy flux in the waves [see e.g. Eq. (2.13a) below], this 
means that the gauge transformation is actually not admissible in the presence of the helical 
symmetry, because the integral of the flux, or total energy emitted, is infinite in that case. 
In the present paper we redo the analysis of the memory term in Eq. (2.2), i.e. 



[6hf] 



memory 



r 



rv rv . 



(2.5) 



and show that it does contribute (together with its iteration at cubic non-linear order) to 
some new log-terms at 4PN and 5PN orders with respect to the results of Paper II. To this 



We pose ft."^ = \/^-g g"^ — i]"^, where g"^ denotes the inverse of the usual covariant metric gaj3, of 
determinant g — det((7Q^), and the Minkowski metric. We use harmonic coordinates throughout, 

i.e. dfsh"^ = 0. 
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end we decompose o"(n, u) into symmetric and trace-free (STF) spherical harmonics: 



+00 



(2.6) 



where (Jl(m) denotes the STF coefficient of order i and ul = Ui^ - ■ ■ rti^. In the Appendix of 
Paper II, it is shown that the log-terms issued from Eq. (2.5) are given by^ 



1=0 



(2.7) 



where di = di^ - ■ ■ and the superscript (— ^ — 3) refers to time antz-derivatives. We 
employ the following convenient notation for a monopolar anti-symmetric wave built from 
an arbitrary function F{u) [25j : 



F{t-r/c) - F{t + r/c) 
2r 



(2.8) 



It was then shown in Paper II that the log-terms (2.7) can formally be rewritten in the form 
of an infinitesimal gauge transformation, i.e. Sh^^'= 2d^°'^^^ — r]°'^dy^"', where the gauge 
vector is explicitly given by 



1 /r 
^« = - In ( - 
^ 2 VA 



El 

£=0 



-cY^'d-d,{a[-'-'^} . 



(2.9) 



"Maximal allowed" gauge tranformation 



However it is crucial to control the occurence of time anti-derivatives in this gauge vector. 
Recall from Paper II that we are looking for so-called near-zone logarithms, present in the 
near-zone expansion when r — )■ of the metric. So we must here discard time ant i- derivatives 



which appear in the near-zone expansion of the gauge vector (2.9). The near- zone expansion 



when r — )■ of the antisymmetric wave (2.8) is regular, and given by 



(2.10) 



Inspection of Eq. (2.9) with the help of the expansion formula (2.10) then shows that the 
monopole i = and dipole i = 1 contributions in (2.9) involve time anti-derivatives. These 



are clearly incompatible with the imposition of the helical Killing symmetry, since for in- 
stance the anti-derivative of the energy flux is the total radiated energy, which is infinite. 
We shall therefore have to define a gauge transformation 77" differing from Eq. (2.9), and 



^ We no longer mention the "memory" origin of this term. For the rest of this Section the notation Sq 
will refer to a modification in the quantity q with respect to Paper II, i.e. which is to be added to the 
corresponding results of Paper II. 
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whose near-zone expansion is free of time anti-derivatives. There is a large number of pos- 
sible choices for such a gauge vector. Here we choose to define 77" by removing from ^'^ a 
minimal number of terms containing the putative anti-derivatives, in such a way that the 
new gauge is still harmonic. This corresponds to a "maximal allowed" gauge transformation 
in the context of helical symmetric spacetimes. Our definition of rj" reads 



° = «°-i{#'-''}''>( 



(2.11a) 
(2.11b) 



where a and (Tj denote the monopolar and dipolar coefficients in the STF multipole expansion 



(2.6), given by (with dfi the solid angle in the direction n) 



ai{u) = 3 I ^nia{n,u) . 



(2.12a) 
(2.12b) 



An easy computation using the definition (2.4) together with the explicit expression (2.3) 
of 2;"^ shows that a{u), which is the average of a{n,u), starts to contribute at order 1/?, 
while the dipole part (5"i(M) starts at order 1/c^. These will correspond to 4PN and 5PN 
terms in the metric, respectively. Here, to be consistent with 5PN, we have to obtain also 
the next-to-leading correction 0{l/c^^) in a, while the leading term is sufficient for &{. We 
find 



cr 



cr,; 



«j ^ ^10 



5c8 

1 



4 

189 



r(4) 



64 
45' 



.(3), 
'ij 



(3) 



+ 



1 

J2 



21 



r{4) 



64 

15^ 



'kl 



+ 



1 



(2.13a) 
(2.13b) 



We recognize the fiuxes of energy and linear momentum carried away by the gravitational 
waves at future null infinity, so that the usual balance equations for the losses of energy E 
and linear momentum P,- read, with such notations. 



dE 
~dt 
d^ 
dt 



12 



cr,: 



(2.14a) 
(2.14b) 



Since the energy fiux is accurate up to next-to-leading order with respect to the quadrupole 



approximation, the quadrupole moment Mij in Eq. (2.13a) has to properly include the IPN 
relativistic corrections. 

Performing the maximal allowed gauge transformation with gauge vector 77" defined by 



Eq. (2.11), we obtain the following extra logarithmic contributions in the quadratic metric 



perturbation (with respect to Eqs. (3.3) in Paper II): 



In 



(2.15a) 
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In 



(2.15b) 
(2.15c) 



It can be checked that dp5h2^ = (modulo some non-logarithmic terms); hence the metric 
is still harmonic. The formula (2.10) then gives the 4PN and 5PN terms in the metric 
perturbation as 



Shi' 



2c2^^ 12c4 



cr 



(2) 



X 

6c^ 



a 



(1) 



In 



IScs""^' 
r 



In 



.14 



+ 

1 



.13 



— 5- (7(2) _ 



j) ^ 18c3 



6iiX a 



In 



A 



O 



1 

34 



(2.16a) 
(2.16b) 
(2.16c) 



As usual this quadratic 4PN + 5PN contribution will generate at cubic order an extra 5PN 
contribution, which is readily computed using techniques similar to those in Paper II. The 
required equation to be integrated is 



A [Shf + Sht] = -^Sh^dijU + -^didhfdiU + O 



1 

J4 



(2.17) 



00 



where U is the Newtonian potential obeying Af/ = —ArcGp^,.^ At 4PN order both 6h' 
and are proportional to a, and thus depend only on time (modulo some non log-terms). 
Using this fact we can readily integrate Eq. (2.17) and obtain the following cubic contribution 
to the log-terms at 5PN order:'' 

2 



6hf + 6hl 



aUln 



A 



+ 



1 

J4 



(2.18) 



Finally, gathering the results we end up with the supplementary contributions {i.e., with 
respect to Eqs. (3.9) and (3.12) in Paper II) of the logarithms at 4PN and 5PN orders in 
the covariant metric as 



. ^tj . 2 - (2) 

— —u o rr cr^ ' 



— -X ' 



6c4 



In 



In 



A 



+ 



.13 



(2.19a) 
(2.19b) 
(2.19c) 



This harmonic metric is rather simple and we have done all computations with it. 



^ Here = miS{x — yi) + rn2S{x — 1/2) is the Newtonian coordinate mass density of the particles. 
^ We also obtain a contribution of far-zone logarithms given by 

However, according to the arguments in Paper II, we can ignore the contribution of far-zone logarithms. 
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2. Cross-checking in a different gauge 



As a check of the computations, we have also used an equivalent metric corresponding to 
a different (non-harmonic) coordinate system, and given by 



'^^00 = 








= - 


G\ 



2G2 



-Ua 



In 



A 



-1 



^aAf/ln 



+ 



.14 



2^ 



13 



(2.20a) 
(2.20b) 
(2.20c) 



This new metric involves a term which is not proportional to ln(r/A) but instead contains 
the logarithm in its "source", originating in a 5PN modification of the source density p^. 
This term, namely 



A- 



aAf/ln ( - 
A 



G^ 



-a 



Grrii 



In 



\x 



yi\ 



\yi\ 

X 



Gm2 

\x - 2/2| 



In 



\y2\ 

X 



(2.21) 



is crucial to take into account in the calculation with the alternative metric ( 2.20[ ). When 
computing the equations of motion and the associated conserved quantities, we must apply 
the metric at the coordinate locations of the particles, so that ln(r/A) becomes ln(|^/i|/A) 
or ln(|7/2|/A), which play the same role as the logarithms generated in Eq. (2.21). In the 



center-of-mass frame, both types of logarithms become in fine In(ri2/A) (modulo irrelevant 
constant terms), where = \yi — y2\ is the coordinate separation. By contrast, in the 
computation with the original metric (2.19), there is no such term a s (|2.21[ ) corresponding 
to a logarithmic modification of the source. Rather, all terms in Eq. ( |2.19 ) are proportional 
to ln(r/A). The two computations are however easily reconciled because the metrics (2.19) 
and (2.20) only differ by a coordinate transformation x'"' = x° + e"(x), namely 



, 2 2 

^fi'oo = ^9oo dteo + — 

c 



dU 



dU 
dy2 



+ 



1 



^9oi = ^9oi - dieo dtei + O 



.13 



where e\ = e^{yA,t), with explicit expression for the gauge vector 



(2.22a) 
(2.22b) 
(2.22c) 



^rV« In 
12c3 VA 



2?^"-'"(a)+'^ 



12 



(2.23a) 
(2.23b) 



The log-terms arising from the non-linear contribution in Eq. (2.22a) are exactly given by 



(2.21), so that this term can be eliminated by the coordinate transformation (2.23), and the 



two computations give the same gauge invariant results. 
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B. Modification of the equations of motion and invariant quantities 



We then compute the modification of the acceleration of one of the particles induced by 
the change (2.19) in the metric. For definiteness we use the harmonic coordinate system 



for this computation. Because the additional logarithmic contributions (2.19) reduce to the 



mere function of time a{t) at 4PN order, the modification of the acceleration occurs only at 
5PN order. We find for the particle 1 (discarding uncontrolled non log-terms) 



5a\ 



-cr ^ ^ n 
ri2C 



6c 



-a. 



(1) 



In 



X J 



(2.24) 



where = {y\ — y^2)/^^2 and v\ = dy\/dt denotes the coordinate velocity. For simplicity we 
no longer write the neglected 6PN remainder. We now check that this modification of the 
acceleration is conservative, in the sense that it corrects the conserved energy E, angular 
momentum Jj, linear momentum Pi, and center-of-mass position Gi of the binary system. 
Indeed, a few calculations (see Paper II for more details) reveal that these corrections, which 
all occur at 5PN order, are given by 



6E 



5Ji 



5Gi 



"2c2 



cr 



q2 



(1) 



6c 



In 



rv2 
X 



r<2 (^2 



In 



ri2 
A 



& 



In 



ri2 
A 



(2.25a) 
(2.25b) 
(2.25c) 
(2.25d) 



where M is the monopole of the source (reducing to m = m\ + m2 at this level of approxima- 
tion). Mi = miy\+m2y\ is the binary's mass dipole, and E is the relativistic 2.5PN-accurate 
binding energy satisfying the balance equation ( |2.14a ). We now work in the center-of-mass 

frame defined by Gi 



0, which gives in particular M, 

'2^ 
G^ 



(1) 



0; hence 



6E 



rriivl + 777,2^2 



SJi = -—aJiln 
G^ 



V A 

fri2 



— o-jM In, , 
6c V A 



SPi 
6Gi = . 

For circular orbits, the extra log-terms in the relative acceleration a* = a\ 



5d 



128 G^ 



m 



n' In 



(2.26a) 

(2.26b) 

(2.26c) 
(2.26d) 
02 reduce to 

(2.27) 



where at this stage we denote r = ri2 and n* = n\2,^ and we have introduced the symmetric 
mass ratio u = mim2/m?, with m = mi + m2 the total mass. The result (2.27) yields 



Note the slight inconsistency in notation here: r — ri2 is the binary's separation in (2.271, while r — \x 



represented the distance to the field point in e.g. Eqs. (2.191 
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in turn a correction in the invariant orbital frequency $7 as a function of the parameter 
7 = Gm/{r(?) (see Paper II for notation) given by 



5^' 



Gm 



'64 



z/^ 7^ In 7 



(2.28) 



where we have used the fact that ln(r/A) = ^\nj (see Paper II). From Eq. (2.28) we obtain 
the correction in 7 as a function of the orbital frequency Q, or rather of the dimensionless 
invariant PN parameter x = [G mVt/ c^Y^^ ^ 



as 



5^ = X 



64 2 5 , 

V X mx 

15 



(2.29) 



It is important to realize that the results (2.27), (2.28) and (2.29) are coordinate dependent. 



They are given here in the harmonic gauge defined by the metric (2.19). Finally, making 



use of these last results, we find the following 5PN logarithmic corrections in the energy and 
angular momentum for circular orbits in the center-of-mass frame:^ 



5E 
5J 



mi' c X 



2 

Gm^ V 



64 



u^x^ Inx 



15 



32 
15 



z/^x^ Inx 



(2.30a) 
(2.30b) 



These last results are coordinate invariant, and we have checked that they come out the 
same with either metric (2.19) or (2.20). The corrections (2.30) are such that 5E = Q6J, 



which is natural because they arise from terms connected to the gravitational-wave fluxes 
[remember Eqs. (2.14)]. These 5PN logarithmic contributions to the energy and angular 



momentum have to be added to the 4PN and 5PN terms already found in Paper II. Note 
that these corrections are proportional to z/^ and therefore affect the results only starting at 
the post-self-force level. 

Next, we compute the redshift observable, which is coordinate invariant for circular orbits 
under the assumption of helical symmetry |22]- In a coordinate system in which the HKV 



reads K^'da = dt + fld^p, it is given by the time component of the four- velocity of the particle: 



Ml 



/9\ -1/2 



(2.31) 



where vf = {c,v\). Inserting the modification of the metric (2.19), together with the stan- 
dard IPN metric for consistent reduction to circular orbits making use of the results (2.28) 



and (2.29), we find the 4PN-I-5PN corrections 



6u\ 



32 



1886 608 



105 105 



A 



1592 
I05" 



z/ X 



x^ In X . 



(2.32) 



Beware that the notation 6E and bJ here has not exactly the same meaning as in Eqs. ( 2.25 )-( 2.26). 
Indeed, there are additional pieces coming from the reduction to circular orbits of the Newtonian parts of 
E and J using Eqs. ( 2.28 1-( 2.29). However, because we have proved in (2.26dl that iGi — 0, there is no 
additional piece coming from the passage to the center-of-mass frame. See also the discussion in Paper II. 
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We introduced the notation A 



[1712 



mi)/m = a/1 — 4i/ for the relative mass difference/^ 



-A. Again, the correction (2.32) occurs 



so that is simply obtained by changing A into 
at the post-self-force level, and does not affect the high-accuracy comparison between the 
post-Newtonian and the perturbative self-force calculations reported in Papers I and II. 
In the following it will be more convenient to work with the inverse of u\, denoted 



13 \ 1/2 



(2.33) 



and which we will still refer to as the redshift observable following Ref . [22] . In terms of this 
variable the correction with respect to paper II reads^^ 



Szi 



32 



2894 80 , 

A 

105 21 



184 
~21 



V X 



(2.34) 



C. Post-Newtonian results for the conserved quantities and redshift observable 



The expressions and J(f2) of the PN binding energy and angular momentum for 

point-particle binaries on quasi-circular orbits have been computed up to 3PN order by 
different groups |26]43T] . More recently, the 3PN expansion of the redshift observable Z\(yC) 
was computed in Paper I, and compared to the numerical result for the gravitational SF in 
linear black hole perturbation theory. 

The dominant logarithmic contributions, which arise at 4PN order [5^ [55] . together 
with the next-to-leading order 5PN contributions, were computed in Paper II. These were 
assumed to only come from gravitational- wave tails, heuristically the scattering of gravi- 
tational radiation by the background curvature generated by the monopole of the source. 
Furthermore, we have now computed some additional 4PN and 5PN logarithmic terms which 
come from the non-linear memory effect, heuristically the gravitational radiation generated 
by the stress-energy distribution of linear waves. All those logarithmic contributions are 
appropriate to conservative helically symmetric PN spacetimes. 



Combining the result (4.12) of Paper II with the new correction (2.30a) obtained above 
we obtain the 5PN-accurate expression for the binding energy: 



E = — mv c X <\ 



+ 



3 

^4 ' 
675 
"6^ 
3969 



12 



X 



11 19 z/2 

1 V \ x" 

8 8 24 



+ 



34445 205 



576 96 

, , 448 

V v tAv] ^ 

128 ^ 15 



-TT 



V — 



45927 



512 



+ z/e5(z/) + 



96 

ulnxj x'^ 
4988 656 



155 2 35 



35 



5184 



ulnx] X 



I x^ 



(2.35) 



where we recall that x = (Gmil/c^)'^^^ (for simplicity we do not indicate the neglected 6PN 
remainder). We introduced some unknown 4PN and 5PN coefficients 64 (z/) and 65 (z/), which 



We assume, without any loss of generality, that mi < m2. 

We take into account the leading- order terms in the expansion as given by u{ 



1+ 
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however are known to be polynomials in the symmetric mass ratio. In Sec. IV B we shall 
be able to obtain precise numerical estimates of 64(0) and 65(0), which encode information 
at leading order beyond the test-particle result (corrections linear in z/). 

The computation of the angular momentum proceeds in the same way. We first derive 
the 4PN and 5PN logarithmic terms due to gravitational- wave tails, similarly to Section IV 



of Paper II, and then add the correction terms given by Eq. (2.30b). The result is 



J 
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1/2 



135 
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+ 
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\ TT^ 

144 24 
64 
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16 
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/ 15309 



I- 



256 



9976 1312 



105 



15 
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(2.36) 



where j/Sy) and j-oiy) denote other unknown coefficients which are also polynomials in v. 
Next, combining the results (4.10) of Paper I and (5.2) of Paper II with the correction term 
(2.32) found above, we obtain the 5PN-accurate expression of the redshift observable m^, or 
rather of its inverse = l/u\, as 
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(2.37) 



The redshift observable Z2 of particle 2 can immediately be deduced from z\ by setting 
A —A. Here P4(z/), q^iy) and "P^iv)-, (I'oiy) denote some still further a priori unknown 
4PN and 5PN polynomials in the symmetric mass ratio. 



Note that in each of Eqs. (2.35)-(2.37) we have added to the usual PN results, valid 



for any mass ratio, the 4PN and 5PN contributions from the test-mass limit for one of the 



The latter point can be proved from the fact that the expression of E (and similarly for J and z\) for 
general orbits, i.e. before restriction to the center-of-mass frame and circular orbits, must be a polynomial 
in the two separate masses m\ and 777-2. 
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particles, known from exact calculations for test particles in the Schwarzschild geometry: 

{ \ — 2x \ 

E = mvc^\ ^ -1 +C'(z/^), (2.38a) 

J = , , , + 0{u') , (2.38b) 

C \/x[l — 6X) 

z^ = y/l-3x + 0{u) . (2.38c) 



D. Post-Newtonian first law of black hole binaries 



In this Section we shall show that the previous post-Newtonian results for the energy 
angular momentum J, and redshift observables zi^2 of binary systems on circular orbits 
satisfy a very important property called the first law of binary black hole mechanics. Here, 
since we are working within the PN framework, the appropriate description of (non-spinning) 
black holes is by structureless point masses; hence we shall also refer to this result as the 
first law of binary point-particle mechanics. 

For this study we shall introduce the total relativistic (ADM) mass of the binary system: 



M = m + 



(2.39) 



where m = mi + m2 is the sum of the post-Newtonian individual masses, i. e. those which 
enter as coefficients of Dirac delta-functions in the stress-energy tensor of point particles; cf. 



Eq. (3.19) below. The various concepts of mass we use will be further specified in Sec. HE 



The ADM mass M, total angular momentum J, and redshifts zi 2, all being given by 



Eqs. (2.35)-(2.37), are functions of the three independant variables of the problem, namely 



the orbital frequency Q that is imposed by the existence of the HKV, and the individual 
masses mi 2 of the particles. We first find, with the above expressions, that the ADM 
quantities obey the usual relation commonly used in PN theory (see e.g. Refs. [26l IM]): 



dM _ dJ 



(2.40) 



This well-known relation is, for instance, extensively used in computations of the binary 
evolution based on a sequence of quasi-equilibrium configurations O |10]. For black hole 
binaries moving on quasi-circular orbits, it states that the gravitational-wave energy and 
angular momentum fiuxes are proportional (with Q being the coefficient of proportionality). 
Here we find that this relation is satisfied for all the terms explicitly computed up to 5PN 
order, including the non trivial 4PN and 5PN logarithmic contributions. 

Largely unrecognized however are the relations which tell how the ADM quantities change 
when the individual masses mi and m2 of the particles vary (keeping the orbital frequency 
Q fixed). By direct partial differentiation of the expressions (2.35)-(2.36) with respect to 
mi and m2, and comparison with (2.37), we find 



dM 



dJ 



drrii drill 



Zl 



(2.41a) 
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From now on we use geometrized units G 
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dM _^dJ 

dmo drrio 



Z2 



(2.41b) 



Again, these relations hold for all the terms we have computed, namely up to 3PN order 



and including the 4PN and 5PN log-terms. Taking further partial derivatives of Eqs. (2.41) 
with respect to the masses yields the simple symmetric relation 



dm2 drrii ' 

which can be viewed as reflecting some "equilibrium" state of the binary system. 



(2.42) 



The three relations (2.40)-(2.41) can be summarized in the first law of the binary black 



hole (or binary point-particle) mechanics 



5M — Vl5J = zi Snii + Z2 61712 



(2.43) 



This first law provides the changes in the ADM mass and angular momentum in response 
to infinitesimal variations of the individual masses of the point particles, weighted by their 
redshift factors. As we shall prove in Sec. [ITTAI that this 1 aw should be correct at all 
FN orders, its verification provides a very powerful test of the intricate post-Newtonian 



calculations yielding Eqs. ( 2.35 )-( 2.37). An interesting by-product of the first law (2.43) is 



the remarkably simple relation 



M — 2QJ = rriiZi + m2^2 



(2.44) 



which can be seen as a "first integral" of the differential law (2.43). It should be stressed 



that the existence of such a simple, linear, algebraic relation between the local quantities 
2:1 2 on one hand, and the global quantities M and J on the other hand, is not trivial. 



In order to prove Eq. (2.44), we take a linear combination of the two equations (2.41) 



make the change of variables (f2, mi, 7712) — ?• (f2, m, z/), where we recall that m 
the total mass and v = 77111712/ m? the symmetric mass ratio, to get 

dM 

m— — = miZi + m2Z2 , 
am 

where we introduced the convenient combination 



mi + m2 is 



(2.45) 



M = M-nj. 



(2.46) 



which can heuristically be viewed as the energy of the binary in a co-rotating frame. Next, 
we notice that the ratio A4/m is dimensionless and symmetric by exchange mi < — > m2 of 
the masses. It must therefore depend on the three independant variables (f2,mi,m2) only 
through the symmetric mass ratio u and the product mQ. This last observation implies that 



d(M/m) B(Mlm) 
m = \l- 



dm 



(2.47) 
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Interestingly, a result equivalent to (2.44) was implicitly derived, in the case of a specific point-particle 



model, at first order in post-Minkowskian gravity: cf. Eqs. (92) and (96) of Ref. |14) . 
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Combining (2.45) with (2.47), and making use of Eq. (2.40), we obtain the relation (2.44). 



Alternatively, we can use Euler's theorem for homogeneous functions: since, from Eq. (2.43) 



the ADM mass M must be a homogeneous function of degree one in J^/^, nii and m2, we 
have^^ 

T X .1/2 dM dM dM , ^ ^ 

M(J,mi,m2) = J^/^TTTTTT^ + mi^^ + maT^, (2.48) 



dmi 



' drrio 



which combined with the first law (2.43) immediately yields the result (2.44). In Sec. IIIC 



we provide a third derivation, based on the integration of the Einstein field equations. 



Note that we have explicitly verified the partial differential equations (2.40) and (2.41) 



and therefore the first law (2.43), only for those terms in the ADM quantities M, J and the 
redshifts Zi 2 which are known, namely for all the terms up to 3PN order and for the log- 
terms occuring at 4PN and 5PN orders. Evidently we could not verify the law for the 4PN 
and 5PN non-logarithmic contributions (and all 6PN and higher order terms) which are left 
out in the PN calculation. However it is not difficult to find the relationships which must be 
satisfied by the polynomials of the symmetric mass ratio parametrizing the unknown terms 



in Eqs. ( 2.35 )-( 2.37). We find that Eq. (2.40) is exactly satisfied at 5PN order if and only 



if the unknown functions ^4(2^), jbii^) in the angular momentum are given in terms of the 
unknowns 64 (z/), 65 (z/) in the energy by 



5 , , 64 
3A{jy) = -^64(2/) + — , 

2 4988 
l<^{iy) = — e'^iv) 



656 
135' 



(2.49a) 
(2.49b) 



Furthermore, we find that the first law is fully satisfied up to 5PN order if and only if the 



unknowns ^4(2^), p^i^) and g4(z/), q^{h') in the redshift (2.37) are also given in terms of the 
functions 64(1'), 65(1/) in the energy, and their i/-derivatives 64(1/), e'r^i^u), by 
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94 (Z^) 
55 (Z^) 



14 
5 



z/e4(z/) 



3969 128 



6 



ue^iiy) + (1 -4z/) q^iu) - 



256 
45927 

1024 



-ly . 



+ 



35 
9976 

945 



ly + 



1312 
13^ 



3 , 3 , , 64 

— ueJu) H 64 z/ , 

28 ' 14 ' 35 ' 

1 , 1 , , 4988 328 

12 ^ 6 ^ ^ 945 45 



(2.50a) 
(2.50b) 
(2.50c) 
(2.50d) 



Thus, all unknown non logarithmic terms at 4PN and 5PN orders are parametrized by only 
two polynomials of the symmetric mass ratio, which can conveniently be chosen to be 64(1^), 



65 (z/) introduced in the binding energy (2.35). Given the high physical significance of the 



first law, as well as the general poof of it given in Sec. Ill (which does not rely on a PN 
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Here the partial derivatives are to be taken using the independent variables (J, 7711,7712) instead of the 
variables (il, 7711,7712) as used in Eqs. ( 2.40 )-( 2.41 ). In this respect notice that 



dM 



dmi 



dM 



J,m2 



dmi 



dJ 



n,m2 



dmi 



£1.7712 



together with a similar expression for 22- 
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expansion), we shall take for granted that it is valid for all PN terms, including the unknown 
ones, so the relationships (2.49)-(2.50) must be exactly satisfied. We will make further use 
of these relations in Sec. lIVl 

In summary, we have proved by a high-order post-Newtonian calculation in this Section 
(building on the results of Papers I and II), that the first law of black hole mechanics (2.43) 
holds in the particular case of a binary system moving on an exact circular orbit. We shall 
further demonstrate, in Sec. Ill, that the law (2.43) is actually a particular case, valid when 
one assumes the existence of a HKV covering the entire spacetime, of the generalized law 
of black hole mechanics obtained by Friedman, Uryii, and Shibata [11]. Before doing so, we 
want to clarify the concepts of ADM and Bondi masses as we use them in post-Newtonian 
theory, i.e. in relation to the relativistic mass M introduced in Eq. (2.39). 



E. Bondi mass versus ADM mass in post-Newtonian theory 

The structure of the gravitational radiation field generated by isolated systems in general 
relativity was elucidated during the early sixties by Bondi, Sachs, Penrose and coworkers, 
who analyzed its asymptotic structure at future null infinity. Of particular importance is the 
Bondi-Sachs mass-loss formula [351 ES], which relates the rate of change of the Bondi mass 
M-q{U) at a retarded instant of time t/,^^ to the gravitational-wave energy flux J^, through 

^ = (2,51) 

where J-" = §y+ |A^p dfi is computed as a surface integral (at future null infinity) of the 
News function M . Another important result was obtained later by Ashtekar and Magnon- 
Ashtekar [37], who showed that (for isolated gravitating systems) the difference between the 
ADM mass Madm and the Bondi mass M-b{U) is equal to the mass carried away by the 
gravitational radiation emitted between the infinite past and the given retarded instant: 

MADM = MB(t/)+ / HU')AU'. (2.52) 

J — oo 



Since the ADM mass is constant, dMAOM/dt/ = 0, the mass-loss formula (2.51 ) immediately 



follows from Eq. (2.52). More generally, similar results hold for the ADM four-momentum 
= (Madm,-P*) and the Bondi four-momentum Pg = (Mb,Pb). 

In addition to the ADM mass and the Bondi mass, many alternative notions of mass 
have been introduced in general relativity (see e.g. [Mj for a recent review). However, while 
the definitions of Madm and Mb only rely on some universal properties of spatial infinity 
and future null infinity, respectively, most of these alternative notions of mass require the 
introduction of one or several additional structure (s) on top of the spacetime metric gap'-, the 
Komar mass Mk, defined in terms of the timelike Killing vector of a stationary spacetime. 



being one specific example (c/. Sec. IIIC). 

In post-Newtonian theory, a background Minkowski metric r^Q,^ is introduced in addition 
to the usual spacetime metric gap- While this, in effect, breaks the manifest general covari- 
ance of the exact theory, the resulting field equations remain covariant under the Poincare 



Here U denotes a null retarded coordinate, which differs from the shorthand u — t — rjc used in Sec. 
where {f, r} are harmonic coordinates. 



II A 
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group of special relativity. In particular, the PN equations of motion of compact binary 
systems (in harmonic coordinates which preserve the Poincare symmetry) are invariant un- 
der time translations, spatial translations, spatial rotations, and Lorentz boosts. Noether's 
theorem then implies the existence of ten conserved quantities associated with these con- 
tinuous symmetries: the post-Newtonian binding energy E, linear momentum P*, angular 
momentum J*, and the vector such that the center-of-mass position reads = K'^ + P'^t. 

The problem of constructing gravitational-wave templates for circularized inspiralling 
compact binary systems in the PN approximation is usually divided into two sub-problems: 
(i) the computation of the center-of-mass binding energy E, using the conservative part of 
the dynamics of the source, and (ii) the calculation of the gravitational-wave energy flux J-", 
associated with the dissipative part of the dynamics, from the wave- zone gravitational field. 
One then postulates the energy balance relation 

^ = -HT) , (2.53) 

which states that the binding energy decreases at a rate determined by the flux of energy 
carried away by the gravitational radiation. The coordinate time T coincides with the proper 
time of an inertial observer far away from the source, where the geometry is essentially fiat. 
In the case of compact binary systems whose components are modelled as point masses, 
both E and J-" have been computed up to 3.5PN order included (see |2] for a review). One 



could in principle verify the energy balance relation (2.53) explicitly from the knowledge of 
all conservative and dissipative contributions to the dynamics of the binary system up to 
3.5PN beyond the leading order radiation-reaction force, which is at leading 2.5PN order; 
this would correspond to 6PN order beyond the Newtonian motion. Since it is, in practice, 
too challenging to compute the PN equations of motion to such high orders, the relation 



(2.53) has to be assumed, rather than derived. Note however that Eq. (2.53) has been proved 



to hold up to the relative 1.5PN order, for a generic matter source IMj ii9j- 



The similarity between the PN energy balance equation (2.53) and the mass-loss formula 



(2.51) is of course striking. However, the mathematical objects entering these two relations 
are conceptually different: for asymptotically fiat spacetimes, the Bondi mass Mb is defined 
in the exact theory as a surface integral at future null infinity, while the binding energy E is 
only defined in the post-Newtonian approximation to general relativity, and its computation 
involves the near-zone PN metric. Although it has not been proved rigorously that the 
two notions of mass coincide, the identification of the Bondi mass with the PN binding 
energy (or, rather, with M = m + E which includes the rest mass m = mi + 7712) seems 
intuitively sound and natural. Henceforth, we shall thus postulate that if U = T — R is 
the asymptotically null outgoing coordinate associated with an asymptotically radiative (or 
Bondi-type) coordinate system {T, R}, then there exists a spacelike hypersurface T = const 
such that 

Mb{U) = M{T) . (2.54) 
The balance equation (2.53) can then be derived from the mass-loss formula (2.51).^^ Con- 



versely, if the PN energy balance is assumed to be valid, then the identification (2.54) must 



hold. The Carter- Penrose diagram depicted in Fig. [T] illustrates the previous discussion. 
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Since Eq. (2.53) is a functional equality, it can also be written as dE/dt = —J-{t), where t is the harmonic- 



coordinate time in the near-zone of the source. 
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FIG. 1: A gravitationally bound isolated matter source is formed at T = Tq, and starts emitting gravita- 
tional radiation. The ADM mass Madm, as computed on that time slice, coincides with the binding energy 
M{To) = m + E{To), which is itself equal to the Bondi mass Mb{Uo), as computed on the asymptotically 
null hypcrsurface U = Uq. At a later time T = Ti, the binding energy M{Ti) = Mq{Ui) has decreased. 
The difference Madm — Mb{Ui) with respect to the constant ADM mass is equal to the energy taken away 
from the source by the gravitational waves emitted between T = Tq and T — Ti, or equivalently between 
J7 = C/o and f7 = Ui. 



While in the exact theory it is not possible to unambiguously split the conservative and 
dissipative parts of the orbital dynamics {e.g. in the binary black hole spacetimes simulated 
in numerical relativity), this can be done in approximation methods such as PN expansions 
and black hole perturbation theory. For example, in post-Newtonian theory, one can discard 
the purely dissipative 2.5PN and 3.5PN radiation-reaction terms, and consider only the 
conservative dynamics at the Newtonian, IPN, 2PN and 3PN orders. Since there is a one- 
to-one correspondance between the near-zone radiation-reaction force affecting the dynamics 
of the source and the fluxes of energy, linear momentum, and angular momentum computed 
from the wave-zone radiation field (see e.g. Refs. |10lllT]), this amounts to considering non- 



radiative spacetimes. From Eq. (2.52) — in which we set J-" = — and the identification 
(2.54), we find that for such PN spacetimes, all the notions of mass that we have considered 
so far coincide, namely 

Madm = Mb = M = const . (2.55) 

See Fig. [2] for an illustration. We shall also assume that a similar result holds for the angular 
momenta, i.e. that the angular momentum J* computed in PN theory from the near-zone 
metric coincides with the total angular momentum of the system, defined at spatial infinity. 



In Sec. II D, the first law of binary point-particle mechanics, as expressed by Eq. (2.43) 



and its first integral (2.44), was derived directly from PN calculations, under the assumptions 
of helical symmetry and asymptotic flatness. When considering only the conservative part 
of the dynamics of a point-particle binary system on a circular orbit, both conditions are 
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FIG. 2: A non- radiative PN spacetime containing an ever-lasting gravitationally bound isolated matter 
source of constant size TZ. The ADM mass Madm coincides with the binding energy M{T) — m. + E{T), 
which is itself equal to the Bondi mass Mb{U), at all times. 



fulfilled by the post-Newtonian metric. For such orbits, the PN results in Sec. II D establish 
algebraically that, for M = m + J and Zi^2 given there, the first law is indeed satisfied 
up to the PN order involved. We shall now show that this first law actually holds more 
generally so that, for the system being considered, it can be expected to be satisfied at all 
(conservative) PN orders. 



III. THE FIRST LAW OF BINARY BLACK HOLE MECHANICS 

Friedman, Uryu and Shibata [H] considered a one-parameter family of solutions of the 
Einstein field equations, describing an arbitrary number of black holes with a generic distri- 
bution of perfect fluid matter sources having compact support. A perfect fluid is entirely 
characterized by its four- velocity m", normalized to gai3u'^u^ = —1, its energy density e{p, s), 
and pressure P{p, s), both functions of the entropy per unit baryon mass s and the conserved 
proper mass density p, such that Vaipu") = 0. 

Each member of the family of geometries {(y'Q,^(A), m"(A), p(A), s(A)} is assumed to have 
a globally defined Killing vector field K". The Noether current associated with K"" assigns 
to each spacetime a conserved charge Q. The main result established in Ref. [H] relates the 
Eulerian variation 6Q = dQ/dA|A=o of the conserved charge Q to the Eulerian variations SAn 
of the horizon surface areas An of the black holes, as well as to the Lagrangian variations 
A(dMb), A(dS') and A(dCa) of the fluid's baryonic mass, entropy, and vorticity. This 



See Ref. [42] for the recent generalization to the magnetohydrodynamical case. 
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generalized first law explicitly reads 



6Q 



[ [/iA(dMb)+f A(dS)+t;"A(dC„)] + (3.1) 



where /€„ are the uniform surface gravities of the black holes, T = T/u^ is the redshifted 
temperature, and = {h — Ts)/u* the redshifted chemical potential (or specific Gibbs free 
energy), with h = {e + P)/ p the specific enthalpy of the fluid. The spacelike velocity field v" 
measures the peculiar velocities of the fluid elements with respect to the Killing vector field 
K"", whose integral curves define some preferred worldlines; it is defined by = M*(i^'° + f"), 
with f^Vot = 0. The matter fields to be varied are given by 

dMb = pM"dS„ , (3.2a) 
d5 = sdMb, (3.2b) 
dCo, = huadM^ , (3.2c) 

with dSfj = —na^/^ d^x the volume element on the spacelike {t = const) hypersurface S, 



covered by Cartesian coordinates {x*}, and over which the integral in Eq. (3.1) is performed; 

is the future-pointing unit timelike vector normal to S, and 7 = det(7jj) the determinant 
of the induced metric 'jij on E. 



A. The first law for point-particle binaries on circular orbits 



We are interested in applying the general result (3.1) to the particular case of a compact 
binary system, whose components move on exactly circular orbits. For such spacetimes, the 
geometry is invariant along the direction of the helical Killing vector K"^ = + f2 where 
n is a constant, which is identified with the angular frequency dip/dt of the circular motion; 
the vectors = (9j)° and 0" = (d^)'^ are part of the natural basis {dt.dr.de.d^), with 
{r, 9, (f)} the spherical coordinate system associated with the Cartesian coordinates {x*} in 
the usual way, i. e. = r sin 9 cos (p, x'^ = r sin 9 sin and x^ = r cos 9. 

For asymptotically fiat spacetimes, the basis vectors t° and 0" are both asymptotic 
Killing vectors, reflecting the invariance of Minkowski spacetime under time translations 
and spatial rotations. Then, the variation 6Q of the conserved charge Q associated with the 
helical Killing vector K'^ is given by [H] 

6Q = 6M -n6J , (3.3) 

where M is the ADM mass, and J the norm of the total angular momentum Jj of the system, 
defined as a surface integral at spatial infinity, in terms of the three rotational Killing vectors 



of the flat metric rjais {cf. Eq. (109) of Ref. IH]). The variational equation (3.3) holds in the 
center-of-mass frame, in which the ADM three-momentum Pj vanishes. 

As discussed in the Introduction, we shall model the compact objects (namely non spin- 
ning black holes and/or neutron stars) as two point particles with "Schwarzschild" masses 
mi and m2. We therefore consider the simple case of a perfect fluid with vanishing pressure 
P, temperature T, specific entropy s, and peculiar velocity f". The redshifted chemical 



potential then simply reads fi = 1/v}. Combining Eq. (3.3) with (3.1), in which we set 
KnSAn = 0, the generalized first law reduces to 

5M-n5J = j 2A(dMb), (3.4) 
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where z = is the redshift factor. Now, in order to evaluate the Lagrangian variation 
A(dMb) of the baryonic mass element (3.2a), we use the explicit expression of the proper 



mass density p for two point masses, namely 



p{x,t) 



I 2 

^ A=l 



(3.5) 



where yA{t) are the coordinate trajectories of the particles {A = 1,2), 6 is the usual three- 
dimensional Dirac distribution, such that J d^x6{x) = 1, and g = det((yfQ,^) is the determi- 
nant of the covariant metric Qais- In a 3 + 1 decomposition, we have y/—g = N^/j, where N 
is the lapse function. Since the four- velocities of the particles are tangent to the HKV 
evaluated at their coordinate locations, i.e. u'X = u\K"{yA,t), we find that (3.2a) explicitly 
reads 



dMb = -d^x J2 



ruA 



A=l 



Na 



5{x - Pa) 



(3.6) 



We have introduced the shorthands Na = N{yA,t) and (Kn) = Kdiff. From the 3-1-1 
decomposition of the time evolution vector in terms of the lapse N and shift N"', namely 
t° = A^n" + A^", the HKV reads = Nn'^ + (A^" + fi0"). Since the shift vector A^" 
and the basis vector 0° are both tangent to S, we get Kafi'^ = —N, which in turn implies 
(Kn)^ = —Na at the location of each particle. This simplifies the expression of the baryonic 



mass element (3.6)^^ which, once substituted into (3.4), yields our final result in the form^^ 



6M — Q6J = zi 6mi + Z2 6m2 



(3.7) 



which we recognize as our first law of binary point-particle mechanics (2.43) derived us- 
ing post-Newtonian theory. This first law compares two neighbouring helically symmetric, 
asymptotically flat solutions of the field equations, and tells how the changes in the "bary- 
onic" masses of the bodies will affect the ADM mass and angular momentum of the binary 
system. Alternatively, the first law of binary point-particle mechanics could probably be 
derived in the Hamiltonian approach by varying the individual degrees of freedom of the 



point masses, namely their positions, momenta and masses. Some consequences of Eq. (3.7) 
have already been explored, notably the existence of the first integral given by 



M — 2QJ = rriiZi + m2Z2 



(3.8) 



In Sec. II D we derived the first law (3.7) directly from post-Newtonian calculations. Here we 



recovered this important relation based on first principles in general relativity; in particular 
our derivation does not rely on a post-Newtonian expansion. Hence we expect the first law 



(3.7), as well as all of its consequences, including Eq. (3.8), to hold at all conservative orders 



in a PN expansion. 



Namely, dMb = p* d^x where = '^a ''^a ^{x — ua) denotes the baryonic coordinate density. 

For clarity and uniformity in notation, we have made the substitution Am^ — > Sm^ for the Lagrangian 

variations of the individual masses mi and m2. 
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B. Analogies with single and binary black hole mechanics 



In this Section we point out some interesting analogies between the result (3.7) and its 



consequence (3.8) for point-particle binaries on one hand, and some well-known relations 
regarding the mechanics (or thermodynamics) of single or binary black holes on the other 
hand. In the particular case of an asymptotically flat vacuum spacetime with two black 



holes on quasi-circular orbits, the general result (3.1) reduces to 
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5M-^15J = Ki — + K2 — . (3.9) 
Stt Stt ^ ' 

This variational relation can be viewed as a generalization to the binary black hole case of 
the celebrated first law of black hole mechanics 5M — VLy{5J = kSA/ISh), which holds for 
any non-singular, asymptotically flat perturbation of a stationary and axisymmetric black 
hole of mass M, intrinsic angular momentum (or spin) J = Ma, surface area A, uniform 
surface gravity k, and angular frequency f^H on the horizon ^Hl HI]. In the binary black 
hole case, the horizon angular velocity (^h of a single rotating black hole is replaced by the 
orbital frequency fl of the binary. 

The surface area A of a black hole is related to its irreducible (or Christodoulou) mass 
mirr through m?j. = A/{167i) |l5lll6]. We may thus substitute kSA/{87i) — )■ (4mi„ k) Smirr 



in both terms in the RHS of Eq. (3.9). Then, comparing the first law of binary black hole 
mechanics (3.9) to Eq. ( |3.7 ) for the analogous binary point-particle case, we notice the 
formal analogies 

m i — > mirr 5 z i — )■ 4mirrK . (3.10) 

Both analogies are rather intuitive and physically appealing. Indeed, one might expect the 
irreducible mass mirr of a non-rotating black hole to be analogous to the "post-Newtonian" 
mass m of a point particle. Furthermore, the surface gravity /t of a black hole may naturally 
be related to the redshift z of the "associated" point mass, via the gravitational redshift, 
or Einstein effect (c/. Detweiler's Gedanken experiment in Sec. II C of |22])- In particular, 
for two Schwarzschild black holes on quasi-circular orbits, but far enough apart so that they 
can be viewed — in first approximation — as isolated, the surface gravity k — )■ (4mirr)~^, 
while for the point particles we have z — y 1. Therefore, in that limiting case the analogies 



(3.10) seem perfectly sound. 



Then, following exactly the same steps as in Sec. IIP but with the irreducible masses of 
the black holes playing the role of the post-Newtonian masses of the point particles, it can 
easily be established that 

M-21]J = Ki — + — . (3.11) 
47r 47r 

This last relation can be viewed as a generalization to the binary black hole case of Smarr's 



formula M - 2nuJ = KA/{A7r) for a single Kerr black hole |17]. The RHS side of Eq. (3.11) 



See Sec. |l] for a discussion of helical symmetry and asymptotic flatness in binary black hole spacetimes 



Notice however that the assumption of helical symmetry, which is necessary to derive the first law (3.9) 



requires that the two black holes are in co-rotation, and must therefore have non-zero spins [51 111). It 



would be interesting to generalize the binary point-particle first law (3.7) by including spin effects, e.g. 
using a pole-dipole model for the spinning point masses. For rotating black holes, the post-Newtonian 
mass m should be identified with the total mass of the black hole, including the effect of the spin S; hence 
= m?„ + S'V(4m2r). See e.g. Ref. [SJ. 
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one for each black hole. This binary 



is the sum of two terms of the form kA/ (An) 
black hole expression should be compared to the RHS of Eq. (3.8) for the point-particle case. 
We then find the analogy 

mz < — > Aml^K , (3.12) 



I.e. 



precisely what is expected from the analogies (3.10), which were drawn from the differ- 



ential relations (3.7) and (3.9). 



It would be interesting to investigate how far the formal analogies (3.10) can be pushed, 



especially in the relevant regime where the two black holes or point masses are in co-rotation 
and interact strongly with each other, namely when 2; 7^ 1. One might for example try to 



recover the point-particle results (3.7) and (3.8) starting from the black hole results (3.9) 



and (3.11), by taking some suitable point-particle limit. In that respect, one might use the 

to derive the gravitational self-force 



method proposed in Ref. [18] to derive the gravitational self-force |19J |50] for the motion 
of a point particle in some background curved spacetime, by scaling down the mass M and 
size R of an extended compact object, while helding its compactness M/R fixed. Such 
an analysis might improve our understanding of the validity of modelling extended objects 
such as black holes by point masses in general relativity. This would be valuable because 
both post-Newtonian methods and gravitational self-force calculations rely heavily on such 
idealizations [21 E] ■ These questions should be addressed in future work. 

Finally, one can also consider the mixed case of a circular-orbit compact binary system 
composed of a black hole and a point particle. This corresponds to the model usually adopted 
in black hole perturbation theory and the gravitational self-force. Calculations similar to 
those detailed above naturally yield the results 



6A 

5M -Vl5J = K V z5m, 

871 

M - 2QJ = K Vmz. 

An 



(3.13a) 
(3.13b) 



C. Alternative derivation of the first integral relations 



In this Section we show that the algebraic first integral relations (3.8) and (3.11) can 
also be obtained by standard techniques derived from the definition of the Komar mass and 
angular momentum (see e.g. [51]). We introduce a Cartesian-type coordinate system {t, x*}, 
as well as the associated spherical coordinate system {t,r,9,(j)}. These coordinates are 
chosen such that the metric gai3 is explicitly asymptotically Minkowskian at spatial infinity, 
i.e. Qap — )■ f]ai3 wheu r — )■ -|-oo. Performing a 3-1-1 decomposition of the four-dimensional 
spacetime, the metric reads 

ds^ = -N^t^ + -fij {dx' + N'dt) {dx^ + N^dt) , (3.14) 

where is the lapse function, A^* the shift vector, and 7jj the three-metric of t = const 
hypersurfaces. Let S be one such spacelike hypersurface, and S = 9S be the two-sphere at 
spatial infinity representing its boundary. 

The notions of Komar mass Mk and Komar angular momentum Jk are commonly intro- 
duced for stationary and axisymmetric spacetimes [52] . Conserved currents can be built from 
the Killing vectors associated with these symmetries, leading to some well-defined notions 
of mass and angular momentum defined as surface integrals over these currents. However, 
the Komar quantities can also be defined for more generic spacetimes, which are neither 
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FIG. 3: A spatial slice E of a circular-orbit compact binary spacetime. The black hole is characterized by 
its horizon area A and uniform surface gravity k, while the point particle has a mass m and redshift factor 
z = and is such that its four- velocity reads u" = u^K". 



stationary, nor axisymmetric, as long as they are asymptotically flat. For such geometries, 
the Komar mass and angular momentum are defined as 



Mk = lim (t> V't'^dSap 



K 



167r 



lim (b V>^d^^^, 



(3.15a) 
(3.15b) 



where the basis vectors = (c?t)" and 0" = (0^)°' are asymptotic Killing vectors, reflecting 
the invariance of Minkowski spacetime under time translations and spatial rotations. The 



integrals in Eqs. (3.15) are performed over a two-sphere Sr of coordinate radius r, which 



goes to the two-sphere at spatial inflnity S as the limit r — )■ -|-oo is taken at the end of the 
calculation.^^ The two-sphere Sr is the boundary of an hypersurface S^, which is itself part of 
the hypersurface S bounded by S. The surface element on Sr reads dSa/s = 2ri^anp]y/a d'^y, 
where n° is the future-pointing unit timelike vector normal to S, and r° the outward- 
pointing unit spacelike vector normal to Sr, and tangent to S. The metric induced on E 
is 7a/3 = gai3 + naUp, and the metric induced on Sr is aajs = 7a/3 — rar/3, with determinant 
cr = det(craf,) in coordinates covering Sr (see Fig. [3]for notations). Note that the integral 



(3.15b) for the angular momentum Jk exhibits the famous "Komar anomalous factor" —1/2 



with respect to the integral (3.15a) for the mass Mk [53] . 

Although and 0" are only asymptotic Killing vectors, our assumption of helical sym- 
metry implies that K"' = -|- 0" is a Killing vector fleld over the entire spacetime. The 
interesting combination of Mk and Jk is therefore given by 



Mk - 2njK = -- 



1 



577 



VK^ dS, 



a/3 



(3.16) 



For stationary and/or axisymmetric spacetimes, and/or are Killing vectors on the entire spacetime 
manifold, and the integrals may as well be performed over any spacelike topological two-sphere enclosing 
the matter source (where the matter stress-energy tensor T"^ ^ 0). 
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Because of Killing's equation V(Q,-ft'^) = 0, the tensor V'^-ft'^ is antisymmetric. Using a 
version of the Stokes theorem for rank 2 antisymmetric tensor fields, the surface integral in 



Eq. (3.16) can be converted into a volume integral over the hypersurface S bounded by 5* 



(and the multiple event horizons T-Ln for spacetimes with black holes): 

/ V'K^ dSap = Y^j ^"^^ + 2 f VpiV'K^) dE« . (3.17) 

The volume element on E reads dS^ = —ria^ d'^x, with 7 = det(7ij) the determinant of the 
spatial metric 7ij. Then, Killing's equation together with the non-commutation of covariant 
derivatives yields the well-known formula V pV^K^ = R^^K^ . Using the Einstein equation 
to replace the Ricci tensor R^p in favor of the stress-energy tensor T^p, we obtain 

Mk - 2nj^ = -^Y.j "^"^^ "^^-P + 2 / - l'^9ap) n'^K^V^d'x , (3.18) 

where T = gajsT""^ . Similar expressions for Mk and Jk (separately) can be found in text- 
books, for stationary and/or axisymmetric spacetimes; see e.g. Refs. [HIIEIIES]- We shall 



now successively apply the general result (3.18) to the particular cases of binary point par- 
ticle and binary black hole spacetimes. 



1. Binary point-particle case 



We consider first a spacetime with two point masses modelling two compact objects on 



a circular orbit; hence we discard the first term in the RHS of Eq. (3.18). The stress-energy 
tensor is generated by two massive worldlines, and reads 



1 ^ 



(3.19) 



A=l 



where yA{t) are the coordinate trajectories of the particles {A = 1,2), with normalized four- 
velocities u% and where za = are the redshift variables. The determinants g = det^ga/s) 
and 7 = det(7jj) are related by y/—g = N Now, remember that the four- velocities 
of the particles are tangent to the helical Killing vector iiT" evaluated at their coordinate 
locations; hence 

uI = u'aK'X, (3.20) 

where we introduced the shorthand K')^ = K°'(yA,t) {cf. Fig. [3]). Inserting (3.19) into (3.18), 
contracting the tensors, and making use of Eq. (3.20), we obtain 



Mk - 2fi Jk 



E 

A=l 



tuaza 



Na 



(3.21) 



where we recall that {Kn) = Kau". As we have already proved below Eq. (3.6) in Sec. Ill A 
above, we have (Kn)^ = —Na at the coordinate locations of the particles. We thus recover 



an expression which is formally identical to our relation (3.8), with the Komar quantities 



playing the role of the ADM mass and angular momentum, namely 



Mk — 2n Jk = rriiZi + m2Z2 



(3.22) 
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Finally, it can be shown that if the foliation (St)^gjg is adapted to axisymmetry at spatial 
infinity (namely n°'(t)a when r +00), then the Komar angular momentum (3.15b) 
coincides with the total angular momentum [HH]: Jk = J ■ On the other hand, for stationary. 



asymptotically flat spacetimes, if the foliation is such that the unit normal coincides with 
the timelike Killing vector at spatial infinity {i.e. N ^ 1 and A^* — when r — i- +00), 
then the Komar mass (3.15a) coincides with the ADM mass [56ll57j. In Ref. [5H], the equality 
Mk = M allowed the definition of a relativistic version of the classical virial theorem for 
stationary, asymptotically flat spacetimes. More recently, Shibata et al. [59] have shown 
that the equality Mk = M holds for a much larger class of spacetimes (in particular they 
could relax the hypothesis of stationarity); see Eqs. (3.7)-(3.9) in Ref. [59] for the required 
asymptotic conditions on the lapse A^, shift A^*, three-metric 7jj, extrinsic curvature Kij, 
and their spatial derivatives. Assuming that our non-radiative, helically symmetric, post- 
Newtonian spacetimes would satisfy these falloff conditions, we then have 



Mk - 2fi Jk = M - 2nj , 



(3.23) 



and we thus recover the algebraic first integral relation (3.8). From this one sees that the 



curious factor of 2 in Eq. (3.8) is actually related to Komar's anomalous factor [53j. 



2. Binary black hole case 



We turn next to the case of a vacuum spacetime with two co-rotating black holes on 
quasi-circular orbits. From the general result (3.18), in which we set T^^ = 0, we find 



Mk - 2fi Jk 



1 

'8^ 



al3 



■Hi 



1 



in 



(3.24) 



H2 



Following Friedman et al. [TT], we write the surface two-form as (\S°'^ = 2K^°'^^'^^/a d?y, 
where is the unique null vector orthogonal to 1-Ln and such that K'^la = —1- For two 
black holes in co-rotation, the helical Killing vector K"' is tangent to the null generators of 
the horizons; hence the surface gravity k is defined by the usual relation K^'S/jsK"' = kK"". 
The integrand in Eq. (3.24) thus reduces to VaKpdS"^^ = —2Ky/a d'^y, yielding 



Mk - 2^] Jk 



1 

An 



KdA + 



Hi 



An 



(3.25) 



H2 



with dA = y/ad^y the surface element. For a single black hole, the zeroth law of black hole 
mechanics states that k is uniform over the event horizon [33] ■ This result was generalized to 
the binary black hole case in Ref. [IT]; hence k can be pulled out of the integrals. Combining 
this with the argument (3.23), we recover the result (3.11). 



On the other hand, several authors have previously established a relation similar to (3.11 ) 
in the case of two co-rotating black holes on quasi-circular orbits, namely [6], |9] 



M -2VLJ 



1 

An 



V„Ard^° 



Hi 



1 

An 



V^NdS"" 



(3.26) 



H2 



where dS" = r" dA. We shall now derive this result from Eqs. (3.23)-(3.24), thus establish- 



ing that for such binary black hole spacetimes, the radial projection of the gradient of the 
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lapse, r^VaN, essentially coincides with the surface gravity k on each horizon, in agreement 
with physical intuition. We assume that the conditions for which (3.23) holds are met. 
Making use of the antisymmetry of the tensor Va-ft'/j, we have 



13' 



-2VaNdS^ 



2KN^np d^° 



(3.27) 



where we performed an integration by parts and used the 3 + 1 decomposition i^" = A^n" + 
m", with m" = N"" + the shift vector in a co-rotating frame. The second term in 
the RHS of Eq. (3.27) can be expressed using the extrinsic curvature tensor Kap, via the 

V/3 In A^. Next, we notice that K^p is both spacelike 



well-known relation Vr 



-K. 



a/3 ■ 



and symmetric, and that ric 



0; hence Eq. (3.24) becomes 



1 



M -2VLJ = ^ V 



, ^ V.AT dS'-^y" ([ m'Ki. dS^ . (3.2^ 



A=l 



Finally, because the HKV field = A^n" + must become null on the event horizons of 
the black holes, we have iV^ = m^m^ on "Hi and 7^2- Following the authors of Refs. |6l |9], 
we impose the usual boundary condition of vanishing lapse on the horizons, and thus find 
that m* = in the second term in the RHS of Eq. ( |3.28 ) (7jj being positive-definite). We 
thus recover the relation (3.26) established in [HI [9], based on the 3+1 decomposition of the 
Einstein field equations, in the presence of a HKV. 



IV. APPLICATIONS OF THE FIRST LAW 



In this Section we consider a number of practical applications of our first law for compact 
binaries. The first is an application when one mass is much smaller than the other, so that a 
perturbative analysis is applicable. The second is an application in which numerical results 



from a perturbation analysis for zi can be combined with relationships like ( 2.49 )-( 2.50) to 
determine previously unknown coefficients in the PN expressions of E and J. 



A. New gauge-invariant energy and angular momentum 

1. Implications of the first law for perturbation analysis 

Perturbation analysis is applicable in the "two-body problem" when the mass of one of 
the bodies is much smaller than the mass of the other. Typically for our purposes, the large 
mass body will be treated as a black hole of mass m2, while the small mass body may be 
treated either as a point particle or as a small black hole of mass mi <^ m2. In this context, 
we have the usual notions of perturbative energy per unit mass and angular momentum per 
unit mass of the particle mi, 

Si = -uit, (4.1a) 
Ji = ui^, (4.1b) 

where uu and Ul^p are covariant components of the particle's four-velocity . It was shown 
in Ref. [22] that zi = l/u\ = Si — QJ'i. As Detweiler pointed out [22], Si and J7i are not 
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separately gauge invariant, while the redshift observable zi is. In the notation of Ref. [22], 
Eqs. (22) and (23) given there are, to first order in the mass ratio q = mi/m2, 



Ji 



r — 2m2 
y^r{r — 3m2) 
r^m2 



y/r - 3m2 



1 - 



r _ _ dh""'' 
r(r — 2m2) 



Arri') 



dr 



(4.2a) 
(4.2b) 



in which use has been made of conditions (A3) and (A4) for quasi-circular orbits, and where 
is a convenience introduced by Detweiler — it does not denote the four-velocity of the 
particle but is defined in (A5) below; and hap is the regularized metric perturbation, which is 
a smooth vacuum solution in the neighborhood of the particle [60].^^ Making use of Eq. (28) 
in [2^ we obtain (see also Eq. (2.7) from Paper I) 



m2 r(l - 3y) _^_odKp 

— = y + — ~a — " ^ 

r b or 



(4.3) 



where we have introduced the dimensionless invariant PN parameter y = (7712^)^/^. It will 
also be convenient to define a coordinate invariant measure of the constant orbit separation, 
via vq = m2/y. Then, using (4.3), Eqs. (4.2) can be rewritten, to first order in q, as 



Jv 



l-2y 

7712 

Vi/(i - Sy) 



1 - 



2 4 



UaU/3 



r(l - 6y) _a_sdh 



dr 12(1 - 2y) 
r(l -2y)_ _ (9/i"^ 



Ay 



dr 



U u' ^ 
or 



12y dr 



(4.4a) 
(4.4b) 



which, as demonstrated at the end of Appendix lAl are not gauge invariant. Nevertheless, 



Eqs. (4.4) can be combined to give 

zi = £i- nji = y/l - 3y (^1 - ^u'^u^h^fi^ 



(4.5) 



where the term in parenthesis is gauge invariant (see [22] and (AlO) in our Appendix [A|). 

Inspired by the present work, new definitions of energy and angular momentum (per unit 
mass) of the particle can be given, which are separately gauge invariant. Indeed, Eq. (2.41a) 
provokes a powerful suggestion to introduce these alternative gauge invariant quantities 



Si 
Ji 



dM 

dm I 
dJ 
dm I 



(4.6a) 
(4.6b) 



In this Subsection, h^p is indeed a regular metric perturbation in the Schwarzschild geometry, and not 



II A 



the PN metric perturbation from flat space used in Sec. [Ill as defined in the first footnote in Sec. 
The partial derivatives with respect to the masses are talccn with constant orbital frequency fi. i'hus we 
really mean for instance £1 = {dM/dmi)\n^rn2- 
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Note that by definition the same combination as in Eq. (4.5) holds, namely 

zi = Si — fijTi . 

Furthermore, the thermodynamic relation for these quantities holds, i.e. 



dJi 
dQ 



(4.7) 



(4.8) 



and (4.8) yields 



This follows from Eq. (|2.40|) and commutation of partial derivatives. Combining Eqs. (4.7) 

dzi 



Ji 



Zi — VL 

dzi 



(4.9a) 
(4.9b) 



The associated invariant energy and angular momentum are defined by Ei = rriiSi and 
Ji = niiJi. Using the first integral relation (2.44), it is straightforward to verify the 



following connections with the ADM mass and angular momentum: 

dM 



El + E2 
J1 + J2 



M + n 
J 



on 
dinj) 



(4.10a) 
(4.10b) 



We now compute the gauge invariant quantities (4.6) or ( 4.9[ ) to first order in the mass 
ratio q. This readily follows from the expressions (4.9) which can be rewritten, using the 
gauge-invariant result for Zi given by (4.5), as (with = m2/y) 

3y) d{u^u^Kp) 



Si 



Ji 



1 - 



u u 



rn(l 
3(1 - 2y) 



drc 



1712 

y(l - 3y) 



1 



3|/ 



drn 



where the final results here expressly hold to first perturbative order, 
show, up to terms O^q"^), that we can write Si = Si + 6S and J7i = ^/i ^ 



6S 



l-2y 



4 



drn 



UaUfS- 



dh''P\ rn{l - Qy) (d{u''u^Kp 



dr 



12 



;i - 2y) V dn 



(4.11a) 

(4.11b) 

These allow us to 
6 J', where 

' dhai3 ^ 



dr 



(4.12) 

together with 6J^ = 6S /Q. Note that the terms in the first brackets have different meanings; 
the first term means the derivative between different circular orbits, and is gauge invariant 
while, as explained in AppendixjA], the second term is not gauge invariant; rather, as pointed 
out in Appendix |A| it means the derivative of spacetime dependence of UaUph^^ at the 
existing quasi-circular orbit with frequency, f2, held fixed. 

It will be interesting to understand the implications of the new gauge invariant quantities 
Si and Ji as we proceed to calculations at second perturbative order. Even without Si and 



Ji, the relations (2.41) have another powerful implication, namely that zi (and hence u\) 



can be calculated directly from the post-Newtonian energy and angular momentum E and 
J coming from the PN equations of motion, instead of the long PN reduction of the defining 



expression (2.31) as was done in Paper I. This thus gives us a new tool through which we 



will be able to compare future PN results that have been obtained in different ways. 
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2. Post-Newtonian expansions of the invariant energy and angular momentum 



The new gauge-invariant energy Ei = rriiSi and angular momentum Ji = miJi of the 
particle 1 (with smallest mass mi) as defined by (4.6) have been obtained at first order 



in the mass ratio q in Eqs. (4.11). This result applies in the strong field regime, and is 



formally valid up to any post-Newtonian order. We now present the PN expressions of £i 
and J7i, valid at a finite PN order but for any mass ratio. These can be obtained using 



Eqs. (4.9) and the PN expansion of zi given by Eq. (2.37). To write PN expressions valid for 
an arbitrary mass ratio, it is appropriate to use the variables x = (mfi)^/^, u = mim2/m^, 
and A = (m2 — mi)/m. We find, up to 3PN order and augmented by the logarithmic 
contributions at 4PN and 5PN orders. 



+ 



+ 



1 

'4 



1 

' 4' 
27 



-A 



27 -. ^ 

\ A + - 

32 32 2 



X + 
19 



3 

16 



16 



-A 



1 



39 2 

Az/ H z/^ + 

16 32 32 



24 
Az/2 



675 675 
256 ^ 256 
"19315 

+ 



A + 
205 



190 205 



9 



-TT 



1728 576 

3969 3969 , 
-A 



-vr 



192 
155 



16 
34445 



72 



X 



205 



X 



1152 



H vr- 

192 



35 



. , o 4865 o 

z/^ H A H 

128 2592 5184 



Az/3 



V 512 

+ 



512 
224 224 



15 



15 



A 



64 
~15 
224 



64^ 7 , , 



7 



Av 

455 
~ 31104^ 

64 
-45' 



-V 



45 



vXyvx 



X 



/45927 45927 



V 2048 



+ 



+ 



+ 



11876 

315 
2494 



2048 
328 



A + z/ 



4988 4988 



Z/2 + 



315 
656 



315 



A~'ip^{v)-'iAq^{y) 



35 



— Az/2 z/" 

15 45 

2494 , 5938 492 , 

H A H V H Av 

35 35 5 



328 



-V 



v\wx\x 



X 



(4.13) 



together with 



1 



+ 
+ 



+ 



1 
2 

27 



2 3 



3 3, 1v 

a; + I - + -A H 

'44 3 



-Av 



27^ 

H A + V ■ 

16 16 

135 135 , 

\ A + 

32 32 



19 



39 



+ 



3863 41 



216 72 
2835 2835 



-vr 



256 256 
64 64 



A 



Az/H z/^ H Av^ - 

16 16 



V 



2 1 2 

-V I X 



X 



304 41 2 

\ vr" 

9 24 



31 
16' 



Az/2 



Z/ + 
973 



6889 41 



\ TT- 

144 24 
7 



v^ + — A z/^ 
324 648 



Av 
91 



3888 



4 1 4 
z/ X 



A + z/ 



64 64 , 

A 

15 15 



10 



10 



-VAsy) 



Aq,i{v) 



64 
— z 
45 



64 



-V 



V \nx ] X 



31 



/ 15309 

n 



+ 



+ 



512 
11876 



A 



15309 
512 



328 . o 
-V + — Az/^ 



4988 4988 , 

+ ^^A-4p5(/^) 



315 15 

9976 9976 

+ 



315 
656 



315 



105 



^ 23752 656 , 

A H z/H Az/ 

105 105 5 



1312 

i 

15 



4Ag5(^^) 



z/ Inx ' 



(4.14) 



The functions Piiy), PbW) and g4(z/), q'5(z/) have been related to the functions e^iy), 65 (z/) 
in the binding energy E by Eqs. (2.50). From these expressions it is easy to verify that 
£1 — VLJi = z\ holds, together with the thermodynamic relation (4.8) and the connections 



to the ADM quantities found in Eqs. (4.10). 



B. Determination of high order PN coefficients in the binding energy 



Black hole perturbation theory is usually formulated as an expansion in powers of the 
mass ratio q = mi/ 1712- However, at first order in g, the symmetric mass ratio v = 5/(1 + 5')^ 
coincides with q, i.e. q = z/ + (9(z/^). In the extreme mass ratio limit g ^ 1, we may therefore 
expand the redshift observable associated with particle 1 as 



zi{x, v) = zschw{x) + uzsf{x) + C(z/^) 



(4.15) 



where we recall that x = (mfi)^/^. The result for a test particle on a circular orbit around 
a Schwarzschild black hole is known in closed form as 2;g|,ji„ = 

VI - 3x; see Eq. (|2^38c|) (for 

simplicity's sake we remove the label 1 indicating the first particle). The invariant relation 
Zsf{x) encoding the first order mass ratio correction was first computed numerically, in the 
Regge- Wheeler gauge, in Ref. [22]. Alternative self-force (SF) calculations based on different 
gauges (Lorenz gauge and radiation gauge) were later found to be in agreement within the 
numerical uncertainties [6T| 162]. 

The conservative gravitational SF effect zsf{x) has also been computed up to high orders 
in the post-Newtonian approximation. From the general form of the near-zone PN metric, 
the conservative SF effect on zi reads 



zsF (x) = ^ 7fc x''^^ + \nx^6, 



(4.16) 



where k is a positive integer, the coefficients 7^ and are pure numbers, the first logarithms 
occur at 4PN order, and the dots stand for terms involving higher powers of logarithms 
[Inxy, with p ^ 2, which are expected not to occur before the very high 7PN order |21j.^^ 

The exact values of the Newtonian, IPN, 2PN and 3PN polynomial coefficients 70, 71, 
72 and 73, as well as those of the 4PN and 5PN logarithmic coefficients ^4 and ^5 can 
immediately be derived from the result (2.37), valid for any mass ratio. These analytical 
results are reported in Table 

Then Paper II showed, for the first time, that it is possible to extract from a SF calculation 
valuable information corresponding to very high orders in the PN approximation (see also 



The general structure of the near-zone PN expansion is known to be of the type ^ a;"/^(lna;)P, where n 
and p are positive integers [53]. For conservative effects n/2 = fc is a positive integer. 
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Papers I and II actually determined the coefficients ak and (3k in the post-Newtonian expansion in powers 
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Coeff. 


Value 


Coeff. 


Value 


ao 


-1 


70 


2 




9 


71 


1 

2 


as 

/35 


-5 

_121 , 41 2 
3 ^32" 
64 
5 

956 
105 


72 
73 
<54 

h 


19 
8 

1621 41 2 
48 32^ 
64 
5 

4988 
105 



TABLE I: The analytically determined PN coefficients {afc,/?^} for ^igp(?/) and {7fc,5fc} for Z'^-p{x). 



Ref. [M] for a similar analysis). Indeed, by fitting the highly accurate SF data for zsF(a;) 
to a PN series of the type (4.16), using the exact values of the analytically determined PN 



coefficients reported in Table [1} the numerical values of the 4PN, 5PN, and 6PN coefficients 
74; 75; 76) and (^e could be determined.^® These are reported in Table [ill Notice in particular 
how the 4PN and 5PN coefficients 74 and 75 could be measured with high precision. Note 
that 74 and 75 coincide with the — )■ limit of the polynomials ^4(1^), (l^iy) and v^iy), Q^i^) 



introduced in Eq. (2.37), namely 



74 = P4(0) + g4(0) + 

75 = P5(0) + ^5(0) + 



1701 

15309 
1024 



(4.17a) 
(4.17b) 



of y = (7712^1)^'^ of tlic self-force effect Ugp on u\ = 1/zi, defined by analogy with. Eqs. (4.15) and (4.16): 
The two sets of coefficients {a^, and {-jk, Sk} can very easily be related using 

X 

zsf{x) = _ - (1 - Sx) ulp{x) . 

For convenience, we provide the values of both sets of coefficients in Tables [l| and [ll[ 
^® The accuracy of the SF data used in Paper II did not allow an unambiguous distinction between the 
effects of the 7PN polynomial {i.e. 77) and logarithmic contributions (67) in the PN expansion of zsf{x). 
However, an 77 (and coefficient is included in Table |ll] since it was used in Paper II to ensure the 
goodness of the fit that was finally obtained for the lower order coefficients. It essentially captures in a 
single term, to the extent available in the data used, the additional contributions from 6j and higher PN 
order coefficients. In the remainder of the present work, we shall disregard any results beyond 6PN order, 
because we believe that a contribution from /?7 confounds a7 |65j . We expect the aj contribution would 
be substantially improved if PN values could be given for the 6PN and 7PN log coefficients. 
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Coeii. 


Value 


Coefi. 


Value 




11/1 7/1 7/ 
— 114. o4 ( 4 / (^Oj 


74 




"5 


-245.53(1) 


75 


-37.72(1) 


Q6 


-695(2) 


76 


+123(2) 


/36 


+339.3(5) 


56 


-311.9(5) 


07 


-5837(16) 


77 


+4210(9) 



TABLE II: The numerically determined values of higher order PN coefficients, based on a fit to 
the SF data given in Paper II. The uncertainties in the last digits are indicated in parenthesis. 



Extending to 6PN order the 5PN-accurate expression of the circular-orbit ADM energy 



M(fi) = m + E{n) given by Eq. (2.35), we have 



M = m 
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(4.18) 



In addition to the 4PN and 5PN unknown coefficients 64(1/) and 65(1/), we introduced further 
6PN unknown coefficients e^i^u) and ejf (i/), which are also polynomials in v. We wish to 
determine the zeroth order coefficients of these four polynomials in the limit — )• 0, i.e. the 
numerical values of 64(0), 65(0), 65(0), and eg'^(O). To do so, we shall use the first law, more 
precisely Eqs. (2.40) and (2.41a), together with the SF results for the redshift observable zi. 
Since Eq. (2.41a) involves the combination M. = M — QJ, it will be convenient to express 



as a function of the energy M alone via the "thermodynamic" relation (2.40). Indeed 
we have d{A4/Q) = Md{l/Q), from which we deduce 



M = 





M{x) 
1^ ' 



(4.19) 



where M is treated as a function of Q or, more precisely, x, holding the individual masses 



nil and m2 fixed. Introducing the PN expansion (4.18) into (4.19) and integrating yields 
the 6PN-accurate expression^^ 
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mv X <1 + 
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Alternatively, we could have subtracted directly the PN expansions (2.35) and (2.36) and use the relations 



(2.49) linking the unknown coefficients ji{v) and i5{v) in the angular momentum to e4^{v) and e^iy). 
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In the test-particle limit, we recover the 6PN expansion of the exact result, which reads 

M = m + miy{y/l - 3x - l) + 0{ 



(4.20) 



(4.21) 



The relation zi = dAi/drrii [see Eq. (2.41a)] establishes a direct link between the redshift 
Zi of particle 1, and the combination Ai = M — QJ of the ADM energy and angular mo- 
mentum. Making the change of variables {Q, mi, 7712) — > (f2,m, z/), this equation becomes 
(remember that mi ^ m2, with A 



m2 



dM 



■ mi)/m 
1 -4z/4 



= 4Z/) 

A dM 



dm 



2m 



(4.22) 



The analogous relation for the redshift Z2 of particle 2 is obtained by changing A into —A. 
We then make another change of variables, namely (fi, m, z/) — > {x,m,v). Since the ratio 



M./m does not depend explicitly on the total mass m but only through x = 
have mdAi/dm = A4 + ^x dA4 / dx; therefore Eq. (4.22) can be rewritten as 



(mn) 



2/3 



we 



m Zl = Ai -\- 



2x dM 
3 dx 



l-Au + AdM 
2 97 



(4.23) 



We can now expand A = y/1 — Au in powers of u, and neglect terms O^u"^) because these 
contributions are not controlled in the self-force calculation of the redshift Zi. Equating the 



terms 0{v) in both sides of Eq. (4.23), we obtain the algebraic relations between the known 



coefficients 7^ and 5k entering the PN expansion (4.16) of the SF effect on zi, and the z/ — )■ 
limit of the unknown coefficients efc(z/) and e)f (z/) entering the PN expansion (4.18) of the 
energy M, which we rearrange as^° 



64(0) 
65(0) 

66(0) 

em 



7 28037 
3 ' 960 ' 


(4.24a) 


9359293 
161280 ' 


(4.24b) 


11 2^ 88209 
3^^+3^^+ 256 ' 


(4.24c) 


(4.24d) 
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Alternatively, the relations (4.24a)-(4.24b) between 64(0), 65(0) and 74, 75 could also be obtained by 



substituting the v ^ Q limit of Eqs. (2.501 into (4.17). 
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Finally, replacing the coefficients 74, 75, 76, and Sq by their known numerical values, as given 
in Table [IT] above, we find 



64(0) = +153.8803(1) , 
65(0) = -55.13(3), 
66(0) = +588(7) , 
(0) = -1144(2) . 



(4.25a) 
(4.25b) 
(4.25c) 
(4.25d) 



Any future post-Newtonian calculation of the 4PN, 5PN, or even 6PN dynamics of point- 
particle binaries will have to be compatible with these numerical values. In principle, 
these could also be recovered from an accurate post-self-force calculation based on second 
order black hole perturbation theory [68] . We now turn to some possible applications of the 



numerical results (4.25) 



C. Other applications left for future work 



We have used the first law (2.43), together with the recently determined numerical values 



of high order FN coefficients in the self-force contribution z^f to the invariant relation Zi{il) 

to compute new FN coefficients in the binding energy E{Q), at 

However, since the first law has been 



[see Eqs. (|4l5|)-(|4l6|)] 

leading order beyond the test-particle approximation 



derived in full general relativity in Sec. Ill A, the relations (2.40)-(2.41) are expected to hold 



at all (conservative) orders in a FN expansion. Making use of these relations, together with 
the numerical results for the exact conservative SF effect 2;sf(^) on the redshift observable 
[22| l6T| [62] . one could in principle compute the fully relativistic SF contribution Esf{^) 
to the binding energy of non-spinning compact binaries [69].^^ This prospect opens up the 
following applications. 



1. ISCO shift induced by the conservative piece of the gravitational self-force 

In the test-particle limit, the innermost stable circular orbit (ISCO) of the Schwarzschild 
geometry is defined as the point of onset of a dynamical instability for circular orbits (the 
circular orbit for which the radial frequency squared of an infinitesimal eccentricity pertur- 
bation turns negative). The orbital frequency of the Schwarzschild ISCO naturally coincides 
with the frequency obtained by minimizing the specific energy -Eschw = (1 ~ 2x) — 3x — 1 



We notice that the numerical values of the coefficients 64(0) and 65(0), as predicted in Ref. using the 
recently computed conservative gravitational self-force correction to the Schwarzschild innermost stable 
circular orbit |67j, are off by ~ f80% and ~ 320%, respectively. This probably reflects the fact that, in 
the extreme mass ratio limit, even a 5PN- accurate formula for the binding energy does not reproduce 
faithfully the exact relativistic result. 
■^^ Note that in order to do so, it must be assumed that the results of a SF calculation, which involve one 
point particle orbiting a Schwarzschild black hole, can be used in conjunction with the first law (2.43), 
which was itself derived for two point masses in the context of PN spacetimes. From this perspective, 
SF calculations would essentially be treated as post-Newtonian calculations formally including all the PN 
corrections (at linear order in the mass ratio). 
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of a test mass in circular orbit around a non-rotating black hole. Going beyond the test- 
particle approximation, the shift of the Schwarzschild IS CO frequency induced by the conser- 
vative part of the gravitational self-force has recently been computed |671[7D]. This genuine 
strong field result has been used extensively as a reference point for comparison with other 
analytical and numerical methods. 

In a post-Newtonian context, one can also compute an ISCO (for any mass ratio) from a 
stability analysis of the conservative part of the PN equations of motion [71] . Alternatively, 
the innermost circular orbit (ICO), or minimum energy circular orbit (MECO), is defined as 
the minimum of the PN binding energy E{x), when it exists [31]. An extensive comparison of 
the SF-induced ISCO shift to the numerous PN-based estimates of ISCO and ICO available 
in the literature was performed in Ref. [66]. In particular, the standard Taylor-expanded 
3PN result [TlJ, based on a stability analysis criterion of the 3PN equations of motion was 
shown to perform astonishingly well. 

Now, it has been shown on very general ground that, for arbitrary mass ratio compact 
binaries, the definitions of ISCO and ICO are formally equivalent [72]; this conclusion does 
not rely on any PN expansion, and only requires that the conservative dynamics of the 
binary system derives from a Hamiltonian.^^ Hence, the exact value of the ISCO frequency 
shift induced by the conservative SF could (in principle) be recovered by minimizing the 
reduced binding energy E/fi = Eschw{x) + u Esf{x) + Oi^u"^) [69]. 

2. Calibration of the potentials entering the effective- one-hody metric 

Within the effective-one-body (EOB) framework, the circular-orbit binding energy E{x) 
is in one-to-one correspondance with the "temporal" potential A = —glf entering the EOB 
effective metric. The knowledge of the self-force correction Eq,y{x) to the test-particle result 
-E'schw(^) thus immediately translates into the knowledge of the coefficient linear in u in the 
potential A,^^ i.e. the function A^^iu) such that A = 1 — 2u + v Aa,^{u) + C(z^^), where 
u = m/r is the usual inverse EOB Schwarzschild-like radial coordinate. 

Furthermore, by combining this result to the recent constraint obtained from a SF/EOB 
comparison of the periastron advance of black hole binaries on circular orbits [23l |6l], one 
could also compute the SF coefficient entering the expression of the "radial" EOB potential 
B = g^^, i.e. the function Bsf{u) such that B'^ = 1 - 2u + f Bsf{u) + 0(i^^). This 
would complete the determination of the two potentials entering the EOB effective metric 
at linear order in the symmetric mass ratio v. Such results would obviously be very useful 
for improving the calibration of EOB models ^1] . 



In a post-Newtonian context however, the location of the ICO needs not agree with that of the ISCO, 
because of the truncation at a finite PN order of the equations defining these notions [TTJ [73] . 
Assuming that the usual mapping between the effective and EOB Hamiltonians holds at all PN orders, 
and that the non-geodesic terms occuring in the expression of the effective Hamiltonian beyond 2PN 
order are proportional to the radial momentum p^, at all PN orders, thus vanishing for circular orbits. 
See Refs. [23l [Tj for more details. 
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3. Comparison with sequences of quasi- circular initial data in numerical relativity 



The resulting expression for the binding energy E{x) could also be used to revisit compar- 
isons with sequences of quasi-circular initial data in numerical relativity. For comparable 
mass black hole binaries, previous comparisons suggested that the convergence of the PN 
series may improve with respect to the extreme mass ratio limit pH [75H77] . This could be 
investigated using the fully relativistic result for E{x), keeping in mind that only the first 
order correction in u beyond the test-particle result would be under control. 

However, higher-order uncontrolled terms with n ^ 2 in E/fi may give only a 

very small contribution to the exact result. Indeed, in a PN expansion, one can check 
that the terms and 0{u^) in Eq. ( |4.18[ ) contribute less than 1% to the total 3PN 

result, up to the Schwarzschild ISCO at x = 1/6. Furthermore, a similar observation has 
recently been put forward in Ref. [77|, for another coordinate invariant relation, namely the 
general relativistic periastron advance A$(x), in the case of non-spinning binary black holes 
on quasi-circular orbits, even in the strong field regime accessible to numerical relativity 
simulations. 
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Appendix A: Details of the self-force analysis 
1. Some relevant definitions 

In Ref. [22], Detweiler defines the quantities Si, TZi, and J7i in terms of the particle's 
four-velocity and four-momentum by 

Si = -uit , TZi^ul, and Ji = ui^ . (Al) 

In general, they are functions of the proper time, s, and the overdot represents d/ds. In 
consequence. 
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Incidentally, we note that the sequence constructed m Ref. [TU] relies on the first law (3.9 1. 
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1 — 2m2/r ' ' ' 

For quasi-circular orbits, Detweiler assumes that £i, IZi and Ji all change slowly, so that 

A ~ 7^l ~ ji ~ 0{q) , (A3) 

while their rates of change vary more slowly still, so that 

Ei^iZi^ Ji^ 0{q^) . (A4) 

Although not entirely necessary, Detweiler found it convenient to define 

being the non-radial part of the four-velocity of the particle, and with all hap terms removed 
{i.e. it corresponds to the non-radial velocity in the background). 

There are two extra quantities that are needed in Sec. IV A and which Detweiler intro- 
duced in his Appendix A of [25]. The first is simply the contraction 

u^u^ha/s = ^ — i-u'^ht,3 + -^u'^h^n , (A6) 

1 — Im^jr 



that follows directly from Eq. (A5), and the second, in which f2 is held fixed, is 

-——u^u'^haf} = --, r^u^htfi irU^Hfi , (AT) 

or [r — 2m2Y r-^ 

in which we depart slightly from [22] , to emphasize that the derivatives do not act on Ex or 
Jx which, as indicated below ( Al ), are functions of the proper time, s, along the orbit; rather 
the derivative is simply acting on components of the metric embedded in the definition ( A5 ). 



2. Gauge transformation properties 



In the application discussed in Sec. |IV A] we are interested in identifying gauge invariants. 
Here we discuss the transformations of m^m^/Iq^, and related quantities, under infinitesimal 
gauge transformations generated by a gauge vector which leads to a metric perturbation 
given by 

A/i,^ = -V,e/3 - . (A8) 

In his analysis [22], Detweiler assumed equatorial symmetry, so that and all its derivatives 
vanish on the equatorial plane, and the resulting metric perturbations are given in his 



Eqs. (B2)-(B7). As indicated in Sec. I A, and emphasized in Sec. II A for quasi-circular 
orbits we require the gauge transformation to obey the implied Killing symmetry, so that 
we must have [dt^Q^d^ = 0. Making use of Detweiler's (B4)-(B6) and the HKV condition, 
and discarding quantities higher than first order, we find 



af5 



[1 — 2m2/rY 



2m') 
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2£i Ji 

- 2m2/r) 



+ 



J? 



2rC + ^d,i 



(A9) 



where the large square brackets contain only metric component transformations. We can 
write this alternatively as 



"Ah 



al3 



+ 



2Si 



1 — 21712/ r \1 



Ji 



2Ji 



\1 — 2m2/r 



21712 /r 
Ji 



1 



(AlO) 



where, in addition, the large round brackets contain terms which vanish to 0{q) at the 
particle. Thus, it is clear that u'^u^/S.hap is gauge invariant when evaluated at the particle. 
This result is needed in (4.5) to establish that zi is gauge invariant through first order in q. 

To evaluate W^u^drAhap at the particle, we need differentiate only the first large round 
bracket, and we obtain the equivalent of Detweiler's equations (B14) and (B15): 



1712 



dr 
dAhap 



dr 



(1 — 2m2/r)2 



C so that 



(1 - Sy) 



(Alia) 
(Allb) 



where the notation |, indicates that the RHS expression has been evaluated at the particle. 
This was used by Detweiler to show that fl [see his Eq. (28)] was indeed gauge invariant. 
To evaluate UaUpdrAh"" , we should differentiate all large round brackets, and we find 
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so that 



(A12b) 



a result not given explicitly in [22j. Notice that neither u"'u"drhap nor UaUpdrh""" is gauge 
invariant, even though u'^u"hai3 = UaUi^h°'" itself is. The terms discussed in (All) and (A12) 
both occur in (4.4). We see that under gauge transformations. Si and J7i transform according 
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to A£i = - 2>y X [d^^t - dt^ci] and AJl = - 2,y x [d^^t - f?t^<^], respectively. 

They are thus radially gauge invariant, but not under Killing-compatible t and (j) changes. 
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